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The influence of a tight isotropic harmonic trap on photoassociation of two ultracold alkali atoms 
forming a homonuclear diatomic is investigated using realistic atomic interaction potentials. Con- 
finement of the initial atom pair due to the trap leads to a uniform strong enhancement of the 
photoassociation rate to most, but also to a strongly suppressed rate for some final states. Thus 
tighter traps do not necessarily enhance the photoassociation rate. A further massive enhancement 
of the rate is found for strong interatomic interaction potentials. The details of this interaction 
play a minor role, except for large repulsive interactions for which a sharp window occurs in the 
photoassociation spectrum as is known from the trap-free case. A comparison with simplified mod- 
els describing the atomic interaction like the pseudopotential approximation shows that they often 
provide reasonable estimates for the trap-induced enhancement of the photoassociation rate even if 
the predicted rates can be completely erroneous. 



I. INTRODUCTION 

Over the past ten years there has been an increasing 
interest in ultracold atomic and molecular physics. This 
interest was stimulated by the successful experimental 
observation of Bose-Einstein condensation (BEG) in di- 
lute atomic gases These atomic condensates exhibit 
many qualitatively new features. Besides their relevance 
to fundamental quantum statistical and possibly even 
solid-state questions a further interesting aspect is that 
the atoms can bind together to form ultracold and even 
Bose-Einstein condensed molecules |^ |^ Q . 

Although so far the only successful way for achieving 
a molecular BEG is based on magnetic Feshbach res- 
onances, alternative schemes are still highly desirable, 
since magnetic Feshbach resonances do not appear to be a 
universal tool. One of the alternative schemes is photoas- 
sociation where two ultracold or Bose-condensed atoms 
absorb a photon and form a bound excited molecule |^|^. 
Although it was demonstrated that this process gener- 
ates cold molecules, the yield is small compared to the 
one obtained by means of magnetic Feshbach resonances. 
The advantage of photoassociation (and related coherent- 
control schemes) compared to Feshbach resonances is, 
however, their assumed wider range of applicability, since 
there is no need for the occurrence of suitable resonances 
and thus no requirement for specific magnetic properties 
of the atoms involved. 

Besides simple one-step photoassociation that yields 
electronically excited molecules there are also resonant or 
non-resonant multi-step schemes leading to the electronic 
and possibly even rovibrational ground state. One of 
the schemes to produce molecules from atoms with help 
of lasers is, e.g., two-color stimulated Raman adiabatic 
passage • The present work discusses only one-photon 
association explicitly, but it is important to note that 
the discussed transition matrix elements are direct in- 
gredients for the modeling of more sophisticated schemes 
like the mentioned two-photon processes. 

Photoassociation is also a powerful tool for the inves- 



tigation of the properties of cold atoms and diatomic 
molecules. The absorption of the photon typically occurs 
at large internuclear distances, and thus the photoasso- 
ciation spectrum provides important information about 
the long-range part of the molecular potential curves as 
well as the collisional properties of atoms 0, 0, ^| . 

The cooling of atomic samples is usually achieved in a 
trap and thus photoassociation experiments in ultra-cold 
atomic gases are performed in the presence of a trap po- 
tential. In most cases these traps are rather shallow so 
that the corresponding trap frequency to is of the order 
of 100 Hz [13. (Here and in the following the frequency 
of the trap corresponds to the one in the harmonic ap- 
proximation). For such a frequency the influence of the 
trap is expected to be negligible. This may, however, 
change for very tight traps. In fact, it was pointed out 
that the atom-molecule conversion process is more effi- 
cient, if photoassociation is performed under tight trap- 
ping conditions as they are, e.g., accessible in optical 
lattices The advantage of using tight confinement 

has stimulated further theoretical investigations and very 
recently some proposals were made that discuss the pos- 
sibility of using the trapping potential itself for the for- 
mation of molecules [ij^ il5|| . 

The study of photoassociation in tight optical lattices 
is of interest by itself, since it is possible to achieve tai- 
lored Mott-Insulator states containing a large number of 
almost identical lattice sites, each filled with exactly two 
atoms The trap frequency of a lattice site in which 
molecules are produced via photoassociation can be of 
the order of 100 kHz The systematic investigation 

of the influence of a tight isotropic harmonic trap on 
the photoassociation process of two alkali atoms form- 
ing a homonuclear molecule is the topic of this work. 
Realistic atom-atom interaction potentials are adopted. 
This allows to check also the range of applicability of 
the (5-function (pseudopotential) approximation for the 
description of the photoassociation process. 

In the pseudopotential approximations the true atom- 
atom interaction is replaced by one that reproduces 
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asymptotically the two-body zero-energy s-wave scat- 
tering. For this choice of the potential ((5-function) 
and if the two atoms are placed in a harmonic trap 
the Schrodinger equation possesses an analytical solu- 
tion [Ulli. The vahdity regime of the pseudopotential 
approximation has been discussed with respect to the en- 
ergy levels for trapped atoms in [2d |. It was shown that 
the use of an energy-dependent pseudopotential (instead 
of the mostly adopted energy-independent one) gives al- 
most correct energy levels for two harmonically trapped 
atoms. Whether this simplified model for the atomic 
interaction is appropriate for the description of photoas- 
sociation in a harmonic trap is, however, not immedi- 
ately evident. Therefore, the present work compares the 
results obtained using realistic atomic interaction po- 
tentials with the ones obtained with either the energy- 
dependent or -independent pseudopotential. 

Photoassociation in tight traps has been studied theo- 
retically before The energy- independent pseudopo- 
tential approximation was adopted and only photoassoci- 
ation into long-range states discussed. Since the present 
work uses realistic atomic interaction potentials, tran- 
sitions to all final vibrational states can be considered. 
This allows to identify two different regimes with respect 
to the influence of a tight trap on the photoassociation 
rate as well as (approximate) rules where a transition 
from one regime to the other is to be expected. 

The outline of this work is the following. First, a brief 
description of the model systems is given in Sec.^ In 
Sec. lIIII thc influence of a tight trap on photoassociation 
is discussed. This includes after a brief general discus- 
sion of the trap influence in Sec. lIII Al the derivation of 
the introduction of an enhance- 



a sum rule in Sec. lIII Bl 

ment or suppression factor in Sec. lIII CI and the discus- 
sion of two regimes in Sec. lIII E)l and IIII El Then the 
case of repulsive atom- atom interactions is considered 
in Sec. lIII F1 The combined influence of trap frequency 
and atom-atom interaction on photoassociation is inves- 
tigated in Sec. lIII CH the validity of the pseudopotential 
approximation in Sec. lIVI Finally, a discussion and out- 
look is given in Sec.0 All equations and quantities in 
this paper are given in atomic units unless otherwise spec- 
ified. 
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II. THE SYSTEM 



Photoassociation of two identical atoms confined in an 
isotropic harmonic trap and interacting through a two- 
body Born-Oppenheimer potential Vint(^) is considered. 
The spherical symmetry and harmonicity of the trap al- 
lows to separate the center-of-mass and the radial inter- 
nal motion ^3 • The eigcnfunctions of the center-of-mass 
motion are the harmonic-oscillator states. Thus the prob- 
lem reduces to solving the Schrodinger equation for the 
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In Eq. J denotes the rotational quantum number, u is 
the harmonic trap frequency, and fi is the reduced mass 
that is equal to m/2 in the present case of particles with 
identical mass m. 

In order to compute the photoassociation spectrum the 
vibrational wave functions ^'(i?)/i? are determined for 
the initial and final molecular states from Eq. (QJ with 
the corresponding Born-Oppenhcimer interaction poten- 
tials ]4nt(^)- The equation is solved numerically using 
an expansion in B splines. For the investigation of the 
influence of the trap on the photoassociation rate Eq. 
is solved for w ^ 0. 

The photoassociation processes most relevant to ex- 
periments on ultracold alkali atoms correspond to tran- 
sitions from two free ground-state atoms (interacting via 
the ground triplet or singlet potential) to the different 
vibrational levels of the flrst excited triplet or singlet 
state 122,123, 13- to hyperflne interaction, two alkali 
atoms can also interact via a coherent admixture of sin- 
glet and triplet states. This work starts by considering 
the photoassociative transition between the two triplet 
states a"^E^ and 1^5]+ for ^Li. A corresponding experi- 
ment is, e. g., reported in The generality of the con- 
clusions drawn from this specific example are then tested 
by considering also other atoms ('^Li and "^^K) or modi- 
fying artificially the interaction strength, as is discussed 
below. 

For the short-range part of the a'^S+ molecular poten- 
tial of Li2 the data in [25j are used, including the van der 
Waals coefficients cited therein. In the case of the I'^S^ 
state data for interatomic distances between R 



4.66 ao 

and R = 7.84 aq are taken from |2^ and are extended 
with ab initio values from |27l | for distances between 
R ~ 3.25 ao and R = 4.50 ao and between R = 8.0 ao and 
R = 30.0 ao. The van der Waals coefficients from 
are used. For a E to S molecular dipole transition the 
selection rule is J = J' ± 1. Assuming ultracold atomic 
gases the atoms interact initially in the J' = state of 
the a^E+ potential. The dipole selection rule leads then 
to transitions to the J = 1 states of With the 

given potential-curve parameters a solution of Eq. in 
the absence of a trap {uj = 0) yields for the fermionic ^Li 
atoms 10 and 100 vibrational boimd states for the a'^I]+ 
(J' = 0) and the l'^S+ (J = 1) states, respectively. In 
the case of the bosonic ""Li atoms there are 11 and 108 
vibrational bound states for the a"^E+ (J' = 0) and the 
I'^Sg (J = 1) states, respectively. 

The electronic dipole moment D{R) for the transition 
a'^Sj 1 '^^ff calculated with a con- 

figuration interaction (CI) method for the two valence 
electrons using the code described in [33| ■ The core elec- 
trons were described with the aid of the Klapish model 
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potential with the parameters given in |31l and polariza- 
tion was considered as discussed in p32| . The resulting 
D{R) (and its value in the separated atom limit) is in 
good agreement with literature data |23, |U 01 • 

In the limit of zero collision energy the interaction be- 
tween two atoms can be characterised by their s-wave 
scattering length age- Its sign determines the type of in- 
teraction (repulsive or attractive) and its absolute value 
the interaction strength. For a given potential curve the 
s-wave (J' = 0) scattering length can be determined us- 
ing the fact that at large distances the scattering wave 
function describing the relative motion (for u; = and 
very small collision energies) reaches an asymptotic be- 
havior of the form 



*s(i?) - sin[fc(ii' 



(2) 



In the present numerical approach discretised continuum 
states arc obtained, since the wave-function calculation 
is performed within a finite R range, i. e. in the inter- 
val [0,i?niax]- Only wave functions that decay before 
or have a node at i?max are obtained. From the anal- 
ysis of the lowest lying discretised continuum state the 
scattering length a^c is obtained by the use of the rela- 
tion ttsc = Rmax — -r with k = -v/z/ZB j35|]. a varia- 
k 

tion of i?max changes the continuum discretization and 
therefore results in i?max-dependent lowest lying contin- 
uum solutions 5* Eo ■ The scattering length extracted from 
converges, however, to a constant value as i?max in- 
creases and Eq approaches zero. Using this method and 
the adopted potential curves the scattering length val- 
ues flsc = — 2030ao and Osc = — 30ao are obtained for 
^Li and ^Li, respectively. These values agree well with 
the experimental ones: Usc = (—2160 ± 250) an (^Li) and 
a^c = (-27.6 ± 0.5) ao (^Li) 

The interaction of two ultracold ®Li atoms is strongly, 
the one of ^Li weakly attractive, as is reflected by the 
large and small but negative scattering lengths. In the 
case of two identical fermionic ^Li atoms the asymmetry 
requirement of the total wave function excludes s-wave 
scattering. Thus the present results are more applicable 
for two ^Li atoms in different hyperfine states (where the 
admixture of a singlet potential would, however, usually 
modify the scattering length), but are actually meant as 
a realistic example for a very large negative scattering 
length, i. e. strong attraction. In order to further check 
the generality of the results also the formation of ^^K2 
is investigated as an example for a small repulsive inter- 
action. In this case photoassociation starting from two 
potassium atoms interacting viathc singlet X^T,'^ ground 
state and transitions into the A^Y,^ state are considered. 
This process is not only experimentally relevant , but 
is at the same time an even further check of the general- 
ity of the conclusions obtained from the investigation of 
the transitions between triplet states in Li2 . 

The data for construc ting the relevant potential curves 
for ■^^K2 are taken from |38ll39ll40j . The resulting poten- 
tial curve for the X^T,'^ state yields a scattering length 



flsc ~ +90 uq. This is in reasonable agreement with the 
experimental value given in |4l] | where Ogc is found to be 
lying between -1-90 og and -f230ao. 

Instead of selecting additional atomic pairs that could 
represent examples for other values of the scattering 
length, the sensitivity of the s-wave interaction on the 
position of the least bound state is used to generate ar- 
tificially a variable interaction strength. The scatter- 
ing length is thus modified by a variation of the par- 
ticle mass. The strong mass dependence of the scat- 
tering length is already evident from its change from 
— 2030ao to — 30ao for the isotopes ^Li and ''Li, respec- 
tively. Experimentally, a strong variation of the inter- 
action strength can be realized by the aid of magnetic 
Feshbach resonances 112, Ii3 ■ 



III. 



PHOTOASSOCIATION IN AN ISOTROPIC 
HARMONIC TRAP 

A. Photoassociation in a trap 



While Eq. ^ yields in the trap-free case (w = 0) both 
bound (vibrational) and continuum (dissociative) states, 
the harmonic trap potential changes the energy spec- 
trum to a purely discrete one, as is sketched in Fig.^ 
Considering the concrete example of two ^Li atoms where 



V(R) 




Internuclear distance R 

FIG. 1: Sketch of the photoassociation process for ®Li2 in the 
presence of a trap (not to scale). With the aid of a laser pho- 
toassociation is induced from the first trap-induced bound state 
{v' = 10) of the a'^Sj state into some vibrational level v of the 
I'^S^ state. Different laser frequencies Ui couple the same initial 
state to diff'crent final states. 

the a'^S]+ state supports the 10 vibrational bound states 
v' = to 9, w' = 10 (J' = 0) denotes the first state 
that results from the trap-induced continuum discretiza- 
tion. This (first trap-induced) state describes the initial 
state of two spin-polarized ^Li atoms interacting via the 
a^^u potential curve, if a sufficiently cold atomic gas 
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in an (adiabatically turned-on) harmonic trap is consid- 
ered. In the present work photoassociation (by means of 
a suitably tuned laser) from this initial state to one of the 
vibrational states v of the potential is investigated 

as a function of the trap frequency lo. In view of the 
already discussed relevant dipole-selection rule the final 
state possesses J = 1 and in the following J' = and 
J = 1 is tacitly assumed. The strength of the photoasso- 
ciation transition to final state v is given by the rate 01 



where X is the laser intensity and 



oc 

J 'f''{R;u;)D{R)^^°'{R;Lj)dR 



(3) 



(4) 



In Eq. g|) 'i'"{R)/R and are the vibrational 

wave functions of the final and initial state, respectively. 
Since the radial density is proportional to it is con- 
venient to discuss 4* instead of the true vibrational wave 
function '^{R)/R. This will be done in the following 
where 5* is for simplicity called vibrational wave function. 
Finally, D{R) is the (i?-dependent) electronic transition 
dipole matrix element between the a'^E+ and the 
state of Li2 introduced in Secl^ D{R) is practically 
constant for R > 25 oq. Eq. (0J is only valid within the 
dipole approximation. The latter is supposed to be appli- 
cable, if the photon wavelength is much larger than the 
extension of the atomic or molecular system. The short- 
est photoassociation laser wavelength corresponds to the 
transition to the highest-lying vibrational state and is 
thus approximately the one of the atomic (2^S — > 2^P 
transition), A = 12680 oq. Although the spatial extent 
of some of the final vibrational states (and of course the 
initial state in the case of shallow traps) has a similar or 
even larger extent, beyond dipole approximation effects 
are neglected in this work. 

The key quantity describing the photoassociation rate 
to different vibrational states v or for variable trap fre- 
quency LO is I "(to) on whose calculation and discussion 
this work concentrates. It is important to note that also 
in the case of more elaborate laser-assisted association 
schemes like stimulated Raman processes that involve 
(virtual) transitions to the v states the transition rate 
is proportional to I^{uj). 

According to Eq. (0J the photoassociation rate depends 
for transitions between long-range states on the Franck- 
Condon factors between the initial and final nuclear wave 
functions, if D{R) is practically constant for large R. In 
the case of alkali atoms the interaction potentials of the 
electronic states can be very long ranged and can con- 
tain numerous rovibrational bound states. Fig. [21 shows, 
e.g., the classical outer turning points i?out of the 100 
(J ~ 0) vibrational bound states of ^Li2 supported by 
the final-state electronic potential curve 1^S+. The or- 
thogonality of the states is achieved by the occurrence of 
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FIG. 2: The classical outer turning points of the vibrational levels 
of the 1 3S+ state of are shown on a linear (solid circles, left 
scale and insert) and on a logarithmic scale (empty circles, right 
scale). 



v' nodes. As v' increases the wavefunctions consist of a 
highly oscillatory short range part with small overall am- 
plitude that covers the range of the v' — 1 wavefunction 
and a large outermost lobe. The I'^S^ state is very long 
ranged, since its leading van der Waals term is —C^/R'^. 
The initial electronic state a'^S+ with leading —Cq/R^ 
van dcr Waals term is shorter ranged. Fig. O shows the 
initial vibrational state for ^Li as a function of the trap 
frequency. This first trap-induced bound state possesses 
V nodes (here v = 10) that arc located in the R range of 
the last trap-free bound state (w = 9). The overall am- 
plitude in this about 25 ao long interval is very small and 
most of the wavefunction is distributed over the harmonic 
trap. 

The squared transition dipole moments P'{u;) are 
shown for ®Li in Fig.^a) for three different trap frequen- 
cies Lu. As mentioned before, the final vibrational levels 
with w > 99 are trap-induced bound states and exist only 
due to the continuum discretization in the presence of a 
trap. If the trap would be turned-off (adiabatically) after 
photoassociation to such a level, the trap induced dimer 
would immediately dissociate (without the need for any 
(radiative or non-radiative) coupling to some dissociative 
state). 

For a fixed trap frequency the photoassociation rate 
generally increases as a function of the final vibrational 
level V, but for small v an oscillatory behavior is visi- 
ble. These oscillations are a consequence of the nodal 
structure of the initial-state wave functions describing 
the atom pair. The 10 nodes (for the shown example of 
^Li) of the initial-state wave function lead to exactly 10 
dips in photoassociation spectrum. Their exact position 
depends on the interference with the nodal structure of 
the final-state wave functions. The oscillatory structure 
of I'"{uj) ends at about u = 55 and beyond that point 
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FIG. 3: Wave functions of the initial a'^Ejt" state of ®Li2 for trap 
frequencies uj = 2it X IkHz (solid), lu = 2it X lOkHz (dashes), uj = 
27r X lOOkHz (chain), and u} = 2n X 500kHz (dots). (The insert 
shows the small R range on an enlarged scale.) 



the rate increases by orders of magnitude, before a sharp 
decrease is observed close to the highest lying vibrational 
bound state {v = 99) . The absence of oscillatory behav- 
ior is a clear signature that for those transitions (in the 
present example for transitions into states with v > 55) 
the Franck-Condon factors are determined by the overlap 
of the outermost lobe of the initial state with the one of 
the final state. 

The comparison of I'"{uj) for the different trap fre- 
quencies shown in Fig.^a) indicates a very systematic 
trend. The transition probabilities to most of the vibra- 
tional bound states increases with increasing trap fre- 
quency. This is in accordance with simple confinement 
arguments, since a tighter trap confines the atoms in the 
initial state to a smaller spatial region. Due to the spe- 
cial properties of harmonic traps, this confinement trans- 
lates directly into a corresponding confinement of the 
pair density (see Eq.([3l)- The probability for atom pairs 
to have the correct separation for the photoassociative 
transition is thus expected to increase for tighter con- 
finements, since a larger Franck-Condon overlap of the 
now more compact initial state with the bound molecu- 
lar final state is expected. However, for the vibrational 
final states close to and above the (trap-free) dissociation 
threshold a completely different behavior is found. In this 
case the photoassociation rate decreases with increasing 
trap frequency, as can be seen especially from the in- 
sert of Fig.^Ja). In fact, a sharp cut-off of the transition 
rate is observed. The transitions to the states that pos- 
sessed the largest photoassociation rate for small trap fre- 
quencies are almost completely suppressed for large trap 
frequencies. Clearly, the simple assumption "a tighter 
trap leads to a higher photoassociation rate due to an in- 
creased spatial confinement" is only partly true. The fact 
that this assumption cannot be valid for all final states 
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FIG. 4: (a) The squared dipole transition moments /"(oj) describ- 
ing transitions from the trap-induced (ii' = 10) initial a^Sj state 
to the vibrational manifold (ti) of the I'^S^ state of '^Li2 are shown 
for the trap frequencies uj = 2it X IkHz (dashes), a; = 27r X lOkHz 
(solid), and uj = 2tt X lOOkHz (dots). The insert shows the transi- 
tions to n = 95 to 100 on an enlarged scale. For a better visibil- 
ity the discrete transitions (marked explicitly with different sym- 
bols in the insert) are plotted as continuous lines, (b) The ratio 
/"(w) [defined in Eq. is shown for lu = 2n X 10 kHz (solid) and 
UJ = 2n X 100 kHz (dashes) as a function of the final vibrational 
level V. (As in (a) the discrete points are connected by lines to 
guide the eye.) The insert shows the transitions to f = 88 to 99 on 
a magnified scale. 



can be substantiated by means of a general sum-rule that 
is derived and discussed in the following subsection. 



B. Sum rule 

Performing a summation (including for lu = an in- 
tegration over the dissociative continuum) over all final 
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vibrational states (using closure) yields 

I{uj) = ^/"(t^) = / *i°'(i?;t^)i^2(i?)*i"'(i?;c^)di?. 

"=o 

(5) 

While the electronic transition dipolc moment D{R) de- 
pends clearly on R for small internuclear separations, it 
reaches its asymptotic value (the sum of the electronic 
dipole transition moments of two separated atoms, Dat) 
at some R value that is much smaller than the typical spa- 
tial extend of the final vibrational states with the largest 
transition amplitudes. (In the example of Li2 this asymp- 
totic limit is reached at about 25 ao.) If the largest pho- 
toassociation amplitudes result from transitions to final 
states whose wave functions are mostly located outside 
this R range, the integral in Eq. |SJl is dominated by the 
R regime in which D{R) is constant. In this case 
can be taken out of the integral and normalization of the 
initial wave function assures /(w) « / = D^^. 

Consequently, for all trap frequencies that are too 
small to confine the atoms into a spatial volume that 
is comparable to the atomic volumes (leading to D{R) ^ 
-Dat) and thus for all traps relevant to this work (and 
presently experimentally achievable) the total dipole 
transition moment / is to a good approximation inde- 
pendent of the trap frequency uj. Therefore, changing the 
trap frequency can only redistribute transition probabili- 
ties between different final vibrational states. Increasing 
the transition rate to one final state must be compen- 
sated by a decrease of the transition probability to one 
or more other vibrational states. 

A conservative estimate of the minimum and maximum 
influence of a harmonic trap on the photoassociation rate 
is obtained from Jmin = -D^i^ and /max = -D^^^, respec- 
tively, where -Dmin (-Dmax) is the minimum (maximum) 
value of the molecular electronic transition dipole mo- 
ment. 

The sum-rule values obtained numerically for the trap 
frequencies shown in Fig.^Ja) are I{uj = 27r x 1 kHz) = 
11.127222, I{uj = 27r X 10 kHz) = 11.12723, and I{uj = 
27r X 100 kHz) = 11.1273. This may be compared to the 
value lim^i^oo \D{R)\'' = Dl, = \D2s-2s + /?2s-2p|' = 
|-D2s-2pP = 11.1272213 obtained from the calculation 
described in Sec. |nj Clearly, the sum-rule lO can also 
be used as a test for the correctness of numerical calcu- 
lations. The very small deviations from the predicted 
sum-rule value may, however, not only be a result of 
an inaccuracy of the present numerical approach, but 
also reflect the (small) R dependence of D{R) that al- 
lows some Lo dependence of the total photoassociation 
rate. This interpretation is supported by the fact that 
the numerically found deviations monotonously increase 
with increasing frequency lo. Larger values of uo lead to a 
spatially more confined [R]^) which in turn probes 
more of the i?-dependent part of D{R). Since D{R) ap- 
proaches D{R = 00) = Dat from above, a small increase 
of / is expected for increasing trap frequencies. As is 



evident from Fig.0](a) (especially the insert), the sum- 
rule fulfillment is achieved by a drastic decrease of the 
photoassociation rate to the highest lying final states. 
This compensates the trap-induced increased rate to the 
lower lying states. Since the rate to the highest lying 
states is orders of magnitude larger than the one to the 
low-lying states, the reduced transition probability of a 
small number of states can easily compensate the sub- 
stantial increase by orders of magnitude observed for the 
large number of low-lying states. 

From Eq. (O it is clear that in those cases where most 
of the contribution to the sum rule stems from the R 
range where D{R) is practically constant, there is also 
no influence of the initial-state wave function. Taking D 
out of the integral yields always the self-overlap of the 
initial-state wave function and thus unity. This is im- 
portant, since it indicates that the sum- rule value is also 
(approximately) independent of the atom-atom interac- 
tion potential. 



C. Enhancement and suppression factor /" 

In order to quantify the effect of a tight harmonic trap 
on the photoassociation rate and to get rid of its variation 
as a function of the final-state vibrational level v (that is 
due to the nodal structure and clearly visible in Fig.^a) 
especially for smaller v) the ratio 



may be introduced. It describes the relative enhancement 
(/^(w) > 1) or suppression {f'"{uj) < 1) of the photoas- 
sociation rate to a specific final state w at a given trap 
frequency to with respect to the reference frequency Wrof ■ 
Although it may appear to be most natural to choose 
the trap-free reference (wicf = 0), a finite value 

offers some advantages. First, a different normalization 
applies to = and uj ^ 0, since in one case it is free-to- 
bound transitions, while it is bound-to-bound transitions 
otherwise. Second, from a numerical point of view it is 
more convenient to treat both cases the same way and 
to avoid the variation of the box boundary i?max that 
would otherwise be necessary for the trap-free case. Fi- 
nally, it may be argued that a non-zero trap reference 
state is in fact more relevant to typical photoassociation 
experiments with ultracold alkali atoms, since most of 
them are anyhow performed in traps. In the present work 
ciJi-cf = 27r X 1 kHz was chosen. This value is sufficiently 
small to represent typical shallow traps in which the in- 
fluence of the trap on photoassociation is supposedly (at 
least to a good approximation) negligible. On the other 
hand, it allows to calculate the transition dipole moments 
with reasonable numerical efforts and thus sufficient ac- 
curacy. 

The ratio /"(w) is shown for two different trap fre- 
quencies LO in Fig.^Jb). For most of the vibrational final 
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states a simple constant regime is observed, i. e. the ra- 
tio /""{uj) is independent of v for all except the highest 
lying states. This constant regime is followed by a rela- 
tively sharp cut-off beyond which the ratio /^'{llj) is very 
small. In the constant regime a 100 kHz trap leads to an 
enhancement by almost 3 orders of magnitude. 

Comparing the results for different lo one notices that 
in the range of final states where a constant behavior 
(with respect to v) is observed, a tighter trap leads to an 
increased photoassociation rate, trap-induced enhanced 
photoassociation (EPA). Due to the cut-off this is, how- 
ever, not true, if the last vibrational states are consid- 
ered. Since the range of constant behavior shrinks with 
increasing trap frequency, there is an increasing range 
of vibrational states for which a tighter trap leads to a 
smaller photoassociation rate compared to a shallower 
trap. In this case trap-induced suppressed photoasso- 
ciation (SPA) occurs. This effect is especially visible 
from the insert of Fig.^fa). The physical origin of the 
two different regimes (constant vs. cut-off) and their lu 
dependence is discussed separately in the following two 
subsections. 



D. Constant regime 

Since in the constant regime the ratio /"(w) is com- 
pletely independent of the final state level v, its value 
(for a given w) and constancy (as a function of v) must 
be a consequence of the influence of the trap on the initial 
state. The initial-state wave function for a ^Li atom pair 
was shown for three different trap frequencies in Fig. |31 
A view on the complete wave function reveals directly 
the confinement of the wave function to a smaller spatial 
volume, if the trap frequency is increased. However, on 
this scale the variation of the wave function at a specific 
value of R appears to be quite complicated. Thus it is 
not at all clear why the enhancement factor has for 
so many states a constant value. A closer look at smaller 
internuclear separations (insert of Fig.l^l reveals that be- 
sides the initial oscillatory part confined to the effective 
range of the atom-atom interaction potential there is a 
relatively large R interval in which the wave functions 
for the different trap frequencies vary linearly with R. In 
this case the slope is very small and the wave function 
is thus almost constant. If the Franck-Condon integral 
is determined only by the value of the initial-state wave 
function in this R window, it produces an almost undis- 
torted image of the final-state wave function. 

However, for the ratio f'"{Lo) this final-state depen- 
dence disappears. The reason is that in the R range 
where the initial-state wave function is almost constant, 
its variation with the trap frequency is also R indepen- 
dent, as can be seen from the insert of Fig. 13 In other 
words, one finds ^'^°'(i?;cj) = C(w) • ^'^°'(i?; w^of)- If no 
final-state dependence occurs, the constant C{lo) is re- 
lated to the ratio / via /c(w) = C^{u>). The validity of 
this argument for the occurrence of a constant ratio / can 



thus be checked (and visualized) in the following way. To- 
gether with the correct wave function ^f^" (i?; uj) the ap- 
proximate one, fc{'-^) • (.R]'^vef)^ is plotted where 
/c(w) is the value of the factor / in the constant regime. 
A convenient way to determine fc{io) follows from the ob- 
servation that the constant regime always starts at v = 0. 
Thus /c(w) = /"^°(w) is the most straightforward way 
of /c(w) determination. In Fig. |Slthe correct wave func- 
tion ^f^^ (R;^) is plotted together with the approximate 
wave function y^^fjui) ■ vj/i" {R;uj^.g() for the trap fre- 
quency uj = 2tt X lOOkHz. The agreement between the 
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FIG. 5: Illustration of the constant photoassociation regime. The 
initial-state wave function iji" (_R;a; = 100) (dots) and an approx- 
imation to it, y/c(a; = 100) ■ <i'^°'(R;uj = 1) (solid), arc shown 
together with the final-state vibrational wave functions for v = 88 
(chain) and v = 94 (dashes). 

two wave functions is clearly good in the shown range 
of R values, but it is better for small R, since at about 
R = 500 ao the two wave functions start to disagree. Be- 
low R = 500 flo the two wave functions agree completely 
with each other, even in the very short R range where 
they possess an oscillatory behaviour. 

The key for understanding the occurrence of the con- 
stant regime is that a variation of the trap frequency 
modifies the spatial extent of the initial-state wave func- 
tion, but leaves its norm and nodal structure preserved. 
As a consequence, the wave function changes qualita- 
tively only in the large R range, while in the short range 
only the amplitude varies (by factor C{lli)). The reason 
for this behaviour is that at small R the wave function is 
practically shielded from the trap potential by the dom- 
inant atom-atom interaction. 

Fig.Elshows also two final-state wave functions {v — 88 
and 94). According to Fig.0|(b) the transition to u = 88 
belongs still to the constant regime {f^^{uj) « fdu;)), 
though to its very end. The transition to w = 94 does 
on the other hand not belong to this regime, since for 
the considered trap frequency f^'^iu;) < fd'^)- As is 
evident from Fig. [S] a constant ratio f "{uj) is observed 



as long as the final-state wave function v is completely 
confined within an R range in which the approximation 

« C{uj) ■ *^"'(i?;wrcf) is well fulfilled. This is 
(for u! = 2tt X lOOkHz) the case for v = 88 for which the 
wave function is confined within R < 600 ao, but not for 
w = 94 whose outermost lobe has its maximum at about 
1350 flQ. Since the R range in which the initial-state wave 
function can be approximated in the here discussed fash- 
ion decreases with increasing trap frequency, the range 
of V values for which f^{uj) ~ fc{^) is valid diminishes 
with increasing trap frequency. 

The following rule of thumb is found to determine those 
vibrational levels v for which the relation f'"{Lo) k, fd^^) 
starts to break down. For trap frequencies and W2 
(with UJ2 > uji) one may define a difference A that quan- 
tifies the deviation of C ■ ^^°'{R;uji) and {R;lo2) 
as A(i?) = C- {R;uji) - ^i"^' (i?;w2) where 

C = V fc{^2)- For example, in Fig. the difference 
A(i?) is the distance between the solid curve and the 
dotted one. The relation f^[uj) « /c(w) breaks down for 
those final states v whose classical turning point lies be- 
yond i?o- -Ro itself is determined by A(i? > Rq) > 10~^. 
In other words, if the last lobe of the final wave function 
overlaps substantially with a region where the deviation 
defined by A is larger than about 10"'^, a clear deviation 
from the constant regime is to be expected. 



E. Cut-off regime 

Once the constant regime of the ratio /" (for a given 
trap frequency) is left, /" is steadily decreasing with v, 
as is apparent from Fig.0|(b). The photoassociation rate 
displays then a rather sharp cut-off behavior (see insert 
of Fig. 21(a))- The most loosely bound vibrational states 
of the final electronic state have in the trap-free case the 
largest rate but possess a very small one in very tight 
traps. For those high-lying states the wave functions have 
a very highly oscillatory behavior for short R values and 
a large lobe close to the classical turning point. This 
outermost lobe determines the Franck-Condon integral, 
if the initial-state wave function is sufficiently smooth in 
this R range. In Fig. the initial-state wave function 
is shown together with the ones for v ^ 96 and 98 (for 
w = 27r X lOOkHz). 

It is evident from Fig. |B1 that for v = 96 the overlap 
of the initial wave function with the last lobe of the fi- 
nal state is very large. In fact, for this trap frequency 
the overlap reaches its maximum for v ~ 96 and 97 
(see Fig. 0](a)), despite the fact that the trap-induced 
relative enhancement factor /"(w) is small (Fig. ^(b)). 
In the case of w = 98 the transition rate is not only 
clearly smaller than for v = 96 ot 97, but it is also much 
smaller than the rate obtained for the same level at much 
lower trap frequencies (10 or IkHz). Clearly, one has 
r=^^{uj) < 1 and thus for uj = 2n x lOOkHz the level 
ti = 98 represents an example for a trap-induced sup- 
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FIG. 6: Illustration of the cut-off regime for a 100 kHz trap. The 
initial-state wave function {R; 100) (solid) is shown together 
with the two final-state wavefunctions for v = 96 (chain) and v = 98 
(dashes) . 



pressed rate (SPA) in contrast to the usually expected 
enhanced photoassociation in a trap (EPA, /^{uj) > 1). 
From Fig. it is clear that the reason for the small tran- 
sition rate to = 98 is due to the fact that the outermost 
lobe of the t; = 98 state lies mostly outside the R range 
in which the initial-state wave function is non-zero. The 
least bound state (in the trap free case), w = 99, possesses 
an even smaller photoassociation rate, since in this case 
the outermost lobe lies practically completely outside the 
non-zero R range of the initial-state wave function. Due 
to the imperfect cancellation of the oscillating contribu- 
tions from the inner lobes, the photoassociation rate for 
V ~ 99 is very small, but non-zero. 

Increasing the trap frequency even more will confine 
the initial-state wave function to a smaller R range and 
thus SPA occurs for smaller v values. The origin of the 
suppression is in fact a quite remarkable feature, since 
from Fig. it is clear that the trap has practically no 
influence on the final states, even if one considers the 
highest-lying ones that have very tiny binding energies. 
This is still true, if the spatial extent of the final state is 
much larger than the one of the trap potential. This may 
be interpreted as a shielding of the trap potential by the 
molecular (atom-atom interaction) potential. The rea- 
son for the different shielding experienced by the initial 
and the final states is not only due to the fact that the 
former lies above the dissociation threshold, since then 
the photoassociation rate should dramatically increase, 
if transitions into the purely trap-induced bound states 
of the final electronic state are considered. This is, how- 
ever, not the case as can be seen for the states w > 99 
in Fig. 21(a) . The different shielding is due to the inher- 
ently different long-range behaviors of the two electronic 
potential curves describing the initial (a'^E+) and the fi- 
nal (l'^E+) state. If one introduces the crossing point 
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Rc of the long range part of the van der Waals potential 
with the one of the inverted harmonic trap, it is defined 

by equating Cn/R" and —^u'^R^ where C„ is the corre- 
sponding leading van der Waals coefficient. At the point 
Rc the trap potential starts to dominate. For example, 
in the case of the trap frequency uj = 2^ x lOkHz one 
finds Rc ~ +825 ao and Rc ~ +f7700ao for a^T.+ and 
f'^S^ of Li2, respectively. 



(uj) for a repulsive interaction 



In order to check the main conclusions of the results 
obtained for ^Li2 also the formation of ■^^K2 is inves- 
tigated. While for ^Li a photoassociation process be- 
tween triplet states was considered, a transition between 
the X^S+ and the states is chosen for ^^K. In 

contrast to the large negative scattering length of two 
^Li atoms interacting via the a'^E+ potential two ^^K 
ground-state atoms interact via a small positive s-wave 
scattering length. The obtained results for the squared 
transition dipole moments are qualitatively very 

similar to the results obtained for ^Li2. This includes 
the existence of a constant regime of f ^{i^) followed by 
a pronounced decrease for the highest-lying vibrational 
states, the cut-off. The rule of thumb for predicting the 
range of v values for which a constant ratio is observed 
does also work in this case. ^^K2 shows thus trap-induced 
suppressed photoassociation for the highest lying states 
with a sharp cut-off in the (w) spectrum very much like 
^Li2. Therefore, the results are not explicitly shown for 
space reasons. 

For a more systematic investigation of the influence of 
the scattering length Ogc and thus the type of interaction 
(sign of flsc) and its strength (absolute value of Osc) the 
mass of the Li atoms is varied. The mass variation al- 
lows for an in principle continuous (though non-physical) 
modification of age from very large positive to negative 
values. With increasing mass an increasing number of 
bound states is supported by the same potential curve. 
Since age is sensitive to the position of the least bound 
state, even a very small mass variation has a very large ef- 
fect, if a formerly unbound state becomes bound. For ex- 
ample, an increase of the mass of ^Li by 0.3% changes Cgc 
from — 2030ao to about -1-850 oq. The (for ^Li unbound) 
11th vibrational state becomes weakly bound. A further 
increase of the mass increases its binding energy until it 
reaches the value for ^Li. It is also possible to modify 
flsc from — 2030ao to -1-850 cq by lowering the mass of 
^Li. A larger mass variation is required (about 18 %) but 
the number of bound states remains unchanged. In this 
case the large positive value of Ogc indicates that the 10th 
bound state is, however, only very weakly bound and a 
further small decrease of the mass will shift it into the 
dissociative continuum. 

In Fig. [71(a) /"(w) is shown for Osc — +850 oq (achieved 
by a 0.3% increase of the mass) and three different trap 
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FIG. 7: As Fig. I^but the scattering length is artificially changed 
to Use = +850 ao (see text for details). The additional insert in the 
right bottom corner of (a) shows the range = 89 to d = 95 on an 
enlarged scale. 



frequencies as an example for a large positive scattering 
length and thus strong repulsive interaction. The overall 
result is again very similar to the one obtained for a large 
negative scattering length. A tighter trap increases the 
transition rate for most of the states, but there is a sharp 
cut-off for large v. The position of this cut-off moves to 
smaller v as the trap frequency is increased. However, for 
a large positive value of Usc an additional feature appears 
in the transition spectrum: a photoassociation window 
visible as a pronounced dip in the /" spectrum for large 
V. For the given choice of age this minimum occurs for 
V = 92. 

The occurence of the dip for age ^ has been pre- 
dicted and explained for the trap- free case in [i^. fisf and 
was experimentally confirmed jig . Fig. |H| shows the last 
lobe of the final-state vibrational wave function '^^^{R) 
together with the initial-state wave function, both for 
u! — 2tt X 100 kHz. The key for understanding the oc- 
currence of the dip for large positive scattering lengths 
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FIG. 8: For ^hi2 (the scattering length in the initial channel is 
modified to a^c = +850 ao) and for a 100 kHz trap the outermost 
lobe of the final-state wave function [v = 92, long dashes, left scale) 
is shown together with the initial-state wave function {v' = 10, 
solid, right scale). Auxiliary horizontal and vertical dotted lines 
are given to assist the eye in identifying regions with a positive or 
negative overlap of the wavef unctions. 



and its absence for negative ones is the change of sign 
of the initial-state wave function as a consequence of the 
repulsive atom-atom interaction. In fact, in the trap-free 
case the position of this node agrees of course with the 
scattering length. As can be seen from Fig. |H| the tight 
trap moves the nodal position to a smaller value, but this 
shift is comparatively small (about 5 %) even in the case 
of a 100 kHz trap. For negative values of age this node 
appears to be absent, since in this case only the extrapo- 
lated wave function intersects the R axis, but this occurs 
at the non-physical interatomic separation Ry, = — Ogc. 
As a result of the sign change occurring for Qsc > the 
overlap of the initial-state wave function with a final state 
for which the mean position of the outermost lobe agrees 
with the nodal position (i?x) vanishes. The probability 
for a perfect agreement of those two positions is of course 
rather unlikely, but as can be seen from Fig. [71(a) and 01 
where also an approximation for I^(uj = 0) was derived, 
the cancellation can be very efficient. 

It should be emphasised that of course also for Osc < 
a number of dips occur as was discussed in the context of 
Fig. El The difference between those dips and the one dis- 
cussed for Qsc ^ is the occurrence of the latter outside 
the molecular regime. While the other dips are a direct 
consequence of the short-range part of the atom-atom in- 
teraction potential and thus confined (for Li2) to v < 55 
corresponding to i? < SOoq, the dip occuring for S> 
can be located outside the molecular regime. This is even 
more apparent from Fig. [51 where the spectra for four 
different positive values of shown together with 

the one for the (physical) value Ogc = — 2030ao (all for 
uj = 2ttx 100 kHz). The values a^c = +2020 ao, +350 ao, 
+115 oq, and +50 ao were obtained by a mass increase of 
~ 0.3%,~ 0.8%,~ 2%, and ~ 6%, respectively. In agree- 
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FIG. 9: (Color online) Dependence of the squared dipole transition 
moments /"(oj) on the scattering length age for transitions from the 
first trap-induced (v' = 11) initial a^Sj state to the vibrational 
manifold (ti) of the l^S^ state of Li2 in a a; = 27r X 100 kHz trap. 
Using masses slightly larger than the one of ^Li a^c = +50 ag 
(blue), asc = +115 ao (red), a^c = +350 ao (brown), and asc = 
+2020 ao (green) were yielded. For comparison, the result with the 
physical mass (age = — 2030ao, black) is also shown. The insert 
shows the transitions to i) = 94 to 99 on an enlarged scale. 



mcnt with the explanation given above, the position of 
the dip moves continuously to larger values of v as the 
scattering length increases, since the position i?x of the 
last node of the initial state lies close to age- Also the 
positions of the other dips depend on asc, but their de- 
pendence is much weaker and involves a much smaller R 
interval. Clearly, the positions of the dips become more 
stable if they occur at smaller v. 

Noteworthy, the positions of the first 10 dips agree 
perfectly for asc = —2030 and +2020 oq. In fact, both 
spectra are on a first glance in almost perfect overall 
agreement, except the occurrence of the additional dip 
for V = 92. According to the discussion of the sum rule 
m Sec lmBl the total sum / should be (approximately) 
independent of the atomic interaction and thus asc. This 
is also confirmed numerically for the present examples. 
The insert of Fig. reveals how the sum-rule is fulfilled. 
The due to the additional dip missing transition probabil- 
ity is compensated by an enhanced rate to the neighbor 
states with larger v. 

In all shown cases with age > there exist 11 bound 
states in contrast to the 10 states of ^Li (asc = —2030 ao). 
As mentioned in the beginning of this section, it is also 
possible to change the sign of age while preserving the 
number of nodes. The corresponding /" spectra (again 
for a; = 27r X 100 kHz) are shown in Fig.^1 The same 
values of asc as in Fig.j^l (+2020 ao, +350ao, +115ao, 
and +50 ao) are now obtained by a decrease of the mass 
by - 18%, - 17.5%, - 16%, and - 13%, respectively. A 
comparison of the two Figs.l^ and 110! demonstrates that 
the position of the outermost dip (for asc S> 0) depends 
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FIG. 10: (Color online) As Fig. |5] but for a variation of a^c by 
means of a reduction of the mass with respect to the one of ^Li. In 
this case, the number of bound states remains unchanged and the 
transition starts from the v' = 10 a^Sj state. 



for a given u solely on Ogc, while the other dips (in the 
molecular regime) differ when changing the total number 
of bound states from 10 to 11. A comparison of the 
results obtained for age = — 2030ao and +2020 ao with 
10 bound states in both cases shows that most of the 
nodes in the molecular regime are shifted with respect to 
each other in such a way that the v range hosting 10 dips 
for flsc = — 2030ao contains 9 dips for age = — 2030ao. 

Turning back to Fig.[3and the question of the influence 
of a tight trap on the photoassociation rate for age S> 
one notices that the position of the additional dip appears 
to be practically independent of w. As was explained 
in the context of Fig.|Sl the reason is that the position 
of the outermost node depends only weakly on lu. For 
the shown example this shift is even for a 100 kHz trap 
small compared to the separation of the outermost lobes 
between neighboring v states. Therefore, the shift is not 
sufficient to move the dip position away from v = 92. 
However, if Usc is, e. g. increased to +2020 ao the crossing 
point i?x shifts in a 100 kHz trap to about 1500 ao and 
changes thus by « 25 %. In this case the dip position 
moves from w = 95 to 94. It is therefore important to 
take the effects of a tight trap into account, if they are 
used for the determination of age usin g ph otoassociation 
spectroscopy the way discussed in [islJl^. 

In order to focus on the effect of the tight trap it is 
again of interest to consider the ratio /^{lo) introduced 
in Sec. nil (Jl For small but positive values of age the ratio 
/" is structurally very similar to the case age = —2030 ao 
shown in Fig.^(b). A uniform constant regime covering 
almost all v states is followed by a sharp cut-off whose po- 
sition shifts to smaller v as uj increases. A similar behav- 
ior is encountered for age ~ +850 ao and lu = 10 kHz as 
shown in Fig.[7|(b). However, for a tighter trap (100 kHz) 
a new feature appears. In this case the relative enhance- 



ment at the dip position {v = 92) is smaller than in the 
constant regime, but larger for the neighbor states. The 
enhancement factor for v ~ 92 is only « 25 % of /c, while 
the one for v = 93 is « 60 % larger than fc- This results 
in a dispersion-like structure in f^. It should be em- 
phasised that this is again remarkably different from the 
other dips in r"{i^) {v < 55) that show the same (con- 
stant) enhancement factor fc as their neighbor states. 



G. Combined influence of trap and atomic 
interaction 

In view of the very important question how the ef- 
ficiency of photoassociation can be improved, Fig.|3 re- 
veals that besides the use of a tight trap a large scattering 
length is also favorable. The photoassociation rate (away 
from the dips) is enhanced by orders of magnitude, if age 
varies from age = +50 ao to age = +2020 ao! In view of 
the already discussed fact that the results for the overall 
spectrum differ for age > and age < only by the 
position of the dips, it is evident that photoassociation 
(or corresponding Raman transitions) are much more ef- 
ficient, if logcl is very large. 

In order to understand the dependence of the FC fac- 
tors of the vibrational final states on the scattering length 
it is instructive to look at the variation of the initial-state 
wave function with Ogc for large R values. This is shown 
in Fig. for a; = 27r x 100 kHz. While a large attrac- 
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FIG. 11: Wave functions of the initial state describing two Li 
atoms in a trap with frequency uj = 2it X 100 kHz for different 
masses that yield the scattering lengths Gbc = — 2030ao (solid), 
f*sc = -1-50 ao (dots), and asc = -1-2020 ao (dashes). The insert 
shows the small _R range on an enlarged scale. 

five interaction (age ^ 0) leads to a very confined wave 
function for the first trap-induced bound state, a large 
repulsive interaction (age 3> 0) does not only result in a 
node (responsible for the photoassociation window dis- 
cussed above), but also to a push of the outermost lobe 
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to larger R values. This push is of course counteracted 
by the confinement of the trap. However, only the high- 
est lying final states probe the very large R range. As is 
apparent from Fig.|21the final states v < 92 probe almost 
completely the range R < 1000 qq. Within this R interval 
the absolute value of the initial-state wave function in- 
creases with the absolute value of Gsc- As a consequence, 
the corresponding FC factors and should increase with 
|asc|. An exception to this is the already discussed oc- 
currence of the photoassociation window (spectral dip) 
that occurs for a positive scattering length, if the posi- 
tion of the node is probed by the final-state wave func- 
tion. Consequently, one expects for the low-lying final 
states (in fact for almost all except the very high-lying 
ones and the ones at the dip position) that an increase 
of Iflscl leads to an increased photoassociation rate. 

An evident question is of course, whether the enhance- 
ments due to the use of tighter traps and tuning of Ugc can 
be used in a constructive fashion? In order to investigate 
this question, one can introduce another enhancement 
factor 

g {uj,asc) = —, r (7) 

with flscrof = Oag (and Lo-cd = 27r x IkHz as before). 
Clearly, a cut through (?^(w, Osc) for constant Cgc is equal 
to /"(w)- A cut for constant uj describes on the other 
hand the relative enhancement of the photoassociation 
rate as a function of a^c- 

The function g^{uj,asc) depends of course on the vi- 
brational state w, but as was discussed before, most of 
the states show a constant enhancement factor /c. Thus 
it is most interestingly to investigate gc{oJ,asc) that is 
defined as the g function for vibrational states for which 
the relation /" = is valid. This excludes the states 
in the cut-off regime and those at or very close to the 
photoassociation window. In Fig. 1121 aAuj. a^n) is shown 
as a function of a^c for different trap frequencies. The 
important finding is that (7c('i',asc) increases as a func- 
tion of Lj and |asc|- In fact, within the shown ranges of 
and Csc the function ,9c (^7 Qsc) raises by 6 orders of mag- 
nitude, if the maximum values of lo (27r x 100 kHz) and 
Osc (±2000 oo) are considered! A more detailed analysis 
shows that the enhancement is almost equally distributed 
among the two parameters, i. e. a factor 10^ stems from 
the variation of lj and about the same factor from vary- 
ing Osc- Thus the enhancement of the photoassociation 
rate due to the two different physical parameters occurs 
practically independently of each other, at least in the 
rather large parameter range considered. It should be 
emphasized that these ranges are realistically achievable 
in present-day experiments. It is interesting to note that 
this finding is not only very encouraging with respect 
to the possible enhancement of photoassociation rates 
and related molecule production schemes, but it shows 
also that the influence of the parameters scattering (osc) 
and characteristic length scale of an isotropic harmonic 
trap (flho = ■\/l/(/^w)) on the photoassociation process 
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FIG. 12: (Color online) Enhancement factor [see Eq. Q] in 
the constant (i)-independent) regime as a function of age for trap 
frequencies u) = 2tt X IkHz (black), uj = 2tt X 10 kHz (red), ui = 
2-K X 30kHz (green), u} = 2tt X 50kHz (blue), w = 2it X 70kHz 
(purple), and ui = 2n X 100 kHz (brown). 

is very different from the one observed for the energy. In 
energy-related discussions (like the one on the validity of 
the pseudopotential approximation in |23) it was found 
that the ratio |asc/aho| determines the behavior. In the 
present case, both parameters and not only their ratio 
are important. 

IV. PSEUDOPOTENTIAL APPROXIMATION 

The bound state of two atoms in a harmonic trap when 
the atom-atom interaction Vint(^) is approximated by 

in d 

a regularized contact potential — asc5^{R)-^R with 

2/i oR 

energy-independent scattering length Osc, was first de- 
rived analytically by Busch at al. [13 . The bound states 
with integer quantum number rit are expressed as 

mSR) = \-^-"'ARe-f''/^T{^v)U{-u, ^,i?2) , (8) 

with R ~ R/a^o and the characteristic length scale Oho 
of the harmonic trap introduced in the end of the pre- 
vious section. ^ is a normalization constant having the 
dimension of the inverse of the square root of volume 
(see below) and v is an effective quantum number for the 

Ea' 3 

relative motional eigcnstate, v ~ — ^ . The energy 

. 2w 4 . 

eigenvalues are given by the roots of the equation 

r(~x/2 + 3/4) ^ J_ 
r(-a:/2 + l/4) V2^' 

where x — E^'^/lo and ^ = flsc/flho- 

The initial-state wave function \E'^'' {R',(^) of two ^Li 
atoms interacting through the a'^E^ potential and the 
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pscudopotcntial wave function 5'^ with the physical 
(trap-free) value of the scattering length = — 2030ao 
are plotted together in Fig. El for the case of a trap fre- 
quency uj = 2t: X lOkHz. As expected, wave function 
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FIG. 13; Wave functions of the first trap-induced bound state {lu = 
2tt X lOkHz) of two ^Li atoms interacting through the full a^Sj 
potential (solid), a pseudopotential with the energy-independent 
(trap-free) scattering length asc = — 2030ao (dashes), and one with 
the energy-dependent value = — 2872ao (dots). The insert 
shows the short R range on an enlarged scale. 

'i'a^^ fails completely for short internuclear separations, 
since it does not reproduce any nodal structure at all. 
In addition, possesses a wrong behavior at i? = 
where it is non-zero. In the long-range part agrees 
better with the correct wave function. There the main 
difference is an evident phase shift between the two func- 
tions. This phase shift is a consequence of the trap and 
vanishes in the absence of the trap — > 0). The physi- 
cal reason for the phase shift is the non-zero ground-state 
energy in a trap (zero-point energy and motion) due to 
the Heisenberg uncertainty principle. As a consequence, 
the scattering of the two atoms in a trap differs from the 
trap-free case even at zero temperature. On the basis 
of an analysis of the energy spectrum of two atoms in 
a harmonic trap it was found that a pseudopotential ap- 
proximation using an energy-dependent scattering length 
a E leads to a highly improved description of two particles 
confined in an isotropic harmonic trap |20ll47l |. 

While the scattering length is defined in the limit 
iJ ^ 0, an energy-dependent scattering length can be 
introduced by extending its original asymptotic defini- 
tion in terms of the phase shift for s-wave scattering 
6oiE) to non-zero collision energies. This yields aE ~ 
—ta.n.So{E)/k with k = ^2\iE. Clearly, the evaluation of 
So{E) requires to solve the complete scattering problem 
and thus also ge can only be obtained from the knowl- 
edge of the solution for the correct atom-atom interaction 
potential. 

The values of were obtained in the following way. 



After a determination of the ground-state energy of two 
®Li atoms from a full calculation (using the realistic in- 
teraction potential), this energy is used in Eq. © to find 
aE (that is inserted in the equation in place of age). More 
details about this so-called self-consistency approach are 
given in [4^ . In this way an energy-dependent scattering 
length aE = — 2872ao is, e.g., found for two ^Li atoms 
in a trap with frequency w = 27r x lOkHz. The result- 
ing wave function is also shown in Fig. together with 
the correct one and the one obtained for Ogc = — 2030ao. 
Clearly, the agreement with the correct wave function is 
very good for large R. For R > 150 oq the wave func- 
tion obtained for ue = — 2872ao is not distinguishable 
from the correct one. Only in the insert of Fig. 1131 that 
shows the wave functions at short internuclear separa- 
tions one sees a deviation. It is caused by the absence of 
any nodal structure and the wrong behavior at i? 
of the pseudopotential wave function. In fact, at short 
distances the introduction of an energy-dependent scat- 
tering length that corrects the phase shift leads to an 
even larger error compared to the use of Ogc- 

The validity of the pseudopotential approximation us- 
ing an energy-dependent scattering length has been dis- 
cussed before. In it was found that applicabil- 
ity of this approximation depends on the ratios Pe/aho 
and |asc/aho| where /Jg = {2^Ce)'^^'^ is the characteris- 
tic length scale of the interaction potential in the case 
of a leading Cq/R^ van der Waals potential. For two 
^Li atoms in a trap with uj = 2tt x 100 kHz that in- 
teract via the a'^E,^ potential those ratios are 0.02 and 
0.59, respectively. These validity criteria are, however, 
based solely on energy arguments. In other words, if 
those ratios are sufficiently smaller than 1, the energy 
obtained by means of Eq. ^ with Osc should agree well 
with the correct one. In the present example of ^Li the 
ratio between the correct first trap-induced energy E^'^ 
and E"*^*^ obtained with the energy-independent pseu- 
dopotential is E^°' /E°^_^ = 0.96 for uj = 2tt x 10 kHz and 
= 0-92 for uj = 2ttx 100 kHz. By construction, 
the energy E^^ agrees of course completely with E 

In Fig. El /^(tj) obtained when using the pseudopo- 
tential approximation with energy-independent scatter- 
ing length is compared to the spectrum obtained for the 
correct atom-atom interaction, both for a trap frequency 
u! = 2tt X 10 kHz. The two results disagree completely 
for V < 60. For higher lying vibrational states {v > 60) 
the agreement is reasonable. (Note, however, the log- 
arithmic scale.) For the highest lying states {v > 95) 
very good agreement is found even on a linear scale (see 
insert of Fig. I14() . Adopting the energy-dependent scat- 
tering length yields quantitative agreement already for 
V > 75, but again a complete disagreement for v < 60. 

The breakdown of the pseudopotential approxima- 
tion (with energy-independent or dependent scattering 
length) for describing photoassociation to the low-lying 
vibrational states is, of course, a direct consequence of 
the wrong short-range behavior of the pseudopotential 
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FIG. 14: Squared photoassociation transition moments I" (u)) for 
®Li in a = 27r X 10 kHz trap calculated with the molecular in- 
teraction potential (solid) or within the energy-independent pseu- 
dopotential approximation (dashes). 



wave functions (Tig. I13|) . From the definition of it 
follows that the pseudopotential approximation fails, if 
the final-state vibrational wave function has a substan- 
tial amplitude in the R range in which the initial-state 
wave function is strongly influenced by the atom-atom 
interaction. An estimate for this R range is (in the 
present case) the already discussed effective-range pa- 
rameter Pe = (2fiCGy^^. Since for large v the final-state 
wave function is dominated by its outermost lobe whose 
position is in turn close to the classical outer turning 
point -Rout, the pseudopotential approximation should 
be valid for i?out > Pe- In the case of ^Li one finds 
= 62.5 cq. According to Fig. [5] the pseudopotential 
approximation should thus be applicable for v > 70. A 
recently performed photoassociation experiment for ^Li 
considered the transition to u = 59 ^3 ■ For this specific 
example the pseudopotential approximation would pre- 
dict a two times smaller rate than the full calculation for 
w = 27r X 10 kHz. 

The validity of the pseudopotential approximation for 
predicting the photoassociation rates to the high lying 
states can also be used to investigate whether the sim- 
ulation of different scattering lengths by mass scaling is 
senseful. This could be questionable, since a change of 
the mass does not only modify the scattering length (by 
moving the position of the least bound state), but also 
the kinetic energy term. Thus it may be argued that the 
discussed influence of the scattering length on the pho- 
toassociation process could partly also be a consequence 
of the modification of the kinetic energy, at least in the 
case of a substantial mass variation as it was required 
for preserving the number of bound states. Within the 
pseudopotential approximation the scattering length is, 
however, a parameter independent of the mass. There- 
fore, in contrast to the case of the full calculation it is 



possible within the pseudopotential approximation to in- 
vestigate the isolated influence of a variation of the scat- 
tering length (keeping the mass fixed). A corresponding 
analysis confirms that mass scaling can in fact be used 
to simulate a modified atom-atom interaction. 

The pseudopotential approximation was used already 
in [23| for an analysis of the change of the photoassoci- 
ation rate due to a scattering-length modification. The 
investigation concentrated, however, on very high lying 
vibrational states close to or even above the trap-free 
dissociation limit. Since for transitions to those states 
the R dependence of the electronic transition dipolc mo- 
ment can safely be ignored, it is sufficient to concen- 
trate on the Franck-Condon (FC) factors. In Fig. El 
the squares of these factors are shown as a function of 
the scattering length for 90 < v < 98 and trap frequency 
Lu = 2Tr X 100 kHz. As in the pseudopotential ap- 
proximation is used for the initial state, but here the 
final-state wave function is obtained by a full numerical 
calculation while an approach based on quantum defect 
theory (QDT) was used in Furthermore, Na2 was 

considered in |2ll | while it is Li2 in the present study. 

For the states 90 < w < 93 shown in Fig. Iisr a) the 
dependence on a ^c in a 100 kHz trap is very similar to 
the one found in |2lj. The rather regular variation with 
flsc is due to the fact that the final-state wave function 
probes the flat part of the initial-state wave function, as 
can be seen in the insert of Fig. I15r a) where the wave 
function for w = 92 is shown together with the initial- 
state wave function for three different values of age- The 
initial-state wave function varies almost linearly with Ogc 
in the Franck-Condon window of the u = 92 final state. 

for 



According to the discussion in Scc. lIII El for a 100 kHz 
trap the states u > 90 belong to the cut-off regime, but 
for w < 93 the enhancement factor /" is still close to its 
value fc in the constant regime (see Figs. 0] and I?)). The 
minima of the FC^ factors for Csc 3> are a consequence 
of the dip discussed in Scc. lIII Fl Since the nodal position 
i?x moves to larger R if age increases, the minimum in the 
FC^ factors moves to a larger value of Osc if v increases. 
While the pseudopotential approximation is capable to 
predict the existence of the dip occuring for age 3> 0, 
its position is not necessarily correctly reproduced in a 
trap. This is due to the fact that the pseudopotential 
overestimates the trap- induced shift of the position of the 
outermost node. For example, if the mass of Li is varied 
such that flsc = +850 flQ is obtained, a 100 kHz trap shifts 
i?x to « -1-810 ao (Fig.lSJ and the dip occurs at w = 92 
(Fig. 01. Using the pseudopotential approximation (with 
flsc = +850 Co) yields on the other hand R^ ~ +580 ao 
and the dip occurs for v ~ 90. This error in the prediction 
of i?x increases with Ogc- 

The final states 94 < w < 98 whose FC^ factors are 
shown in Fig. I15r b') probe on the other hand the non- 
linear part of the initial-state wave function (close to the 
trap boundary). Consequently, the dependence on Ugc 
differs from the one found in |21(. While for 90 < w < 92 
the FC^ factors arc first decreasing and then increasing. 
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FIG. 15: (a) Squared Franck-Condon factors between the final 
vibrational states v = 90 (solid), 91 (dots), 92 (dashes), or 93 
(chain) of the I'^S^ state and the initial-state pseudopotential wave 
functions as a function of the scattering length age ■ The trap fre- 
quency is a; = 27r X 100 kHz. The insert shows the v = 92 final- 
state wave function together with the pseudopotential wavcfunc- 
tions for age = — 2000ao (dots), age = (dashes), asc = +2000 ao 
(chain), (b) As (a), but for v = 9A (solid), 95 (dots), 96 (dashes), 
97 (chain), v=98 (dot-dash-dash). The insert shows the v = 97 
final-state wavefunction and the pseudopotential wave functions for 
asc = — 2000ao (dots), asc = (dashes), asc = -1-2000 ao (chain). 



if Csc varies from — GOOOoo to -1-6000 oq, the ones of 93 < 
V < 96 arc purely decreasing. For v = 97 and 98 the FC^ 
factors are on the other hand increasing with age- 

In view of the fact that the scattering length of a given 
atom pair may be known (for example from some mea- 
surement), but the corresponding atom-atom interaction 
potential is unknown, it is of course interesting to investi- 
gate whether the pseudopotential approximation allows 
to predict the enhancement factor also in the constant 
regime, i.e. whether it correctly reproduces /c(w). This 
would allow for a simple estimate of the effect of a tight 
trap on the photoassociation rate in the constant regime 
that covers most of the spectrum. In order to determine 



fc{uj) it is sufficient to analyze the ratio of the initial-state 
wave function 4"°^^ for the trap frequencies lo and Wref- 
This comparison may be done at any arbitrary internu- 
clear separation i?iin provided it belongs to the linear 
regime. The result is 



(i?lin;Wrof) 



(10) 



A very special and simple choice that guarantees that 
R\in belongs to the linear regime is i?iin = 0. With this 
value of one finds from the analysis of '5°^^ 



where A{lu) 



pseudo^^-j - 



Aiu;) 



1 2 



cf 



— , (11) 

U! 



IS 



the normalization 
dE 



factor fulfilling 



Depending on the level 



of approximation one may use either Ogc or in the 
evaluation of A. An even simpler estimate of the influ- 
ence of a tight harmonic trap on the photoassociation 
rate is obtained, if the atom-atom interaction potential 
is completely ignored in the initial state. The harmonic- 
oscillator eigenfunctions at i? = i?iin = yield 



3/2 



(12) 



In Fig.^]the enhancement factors fd^^) calculated at 
the different levels of approximation are shown as a func- 
tion of the trap frequency uj. The results obtained for ^Li 
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FIG. 16: Dependence of the enhancement factor fc on the trap 
frequency i/ = lj/(2it) for ^Li2 using the molecular (solid), the 
energy-independent (dots) and -dependent (dashes) pseudopoten- 
tial, or the harmonic-oscillator (chain) wave functions. The insert 
shows the same curves but for ''^K2. 

and '^^K are compared to each other. Remind, in the lat- 
ter case the scattering length Ogc = +90 ao has a much 
smaller absolute value than for ^Li (age = — 2030ao). 
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Consequently, one expects the atom-atom interaction to 
be less important. This is confirmed by (the insert of) 
Fig. El The results obtained for fd^^) with the aid of 
the different approximations discussed above are in very 
good agreement with the correct result in the case of 
•^^K. Even the simple harmonic-oscillator model predicts 
the enhancement factor in the constant regime very ac- 
curately. 

It should be emphasized that the correct prediction of 
the enhancement factor by the simplified approximation 
works, although the prediction of the rates is completely 
wrong fFig. ll4|) in this constant regime (small v). In the 
case of a large absolute value of the scattering length 
(like for ^Li), i.e. for a strong atom-atom interaction, 
the frequency dependence of /c(w) predicted by the sim- 
plified models is on the other hand not very accurate. In 
fact, the simple harmonic-oscillator model clearly over- 
estimates the enhancement factor for large lo. The pseu- 
dopotential approximation yields much better results, es- 
pecially if the energy-dependent scattering length is 
used. (As already mentioned, qe is, however, only avail- 
able from the knowledge of the exact atom-atom interac- 
tion.) 

In view of the usefulness of Eq. Hll() for obtaining an 
estimate of the enhancement factor fc{i^) but the rather 

dE 

complicated procedure to calculate -jr— required for ob- 

taining A{uj), it is interesting to test whether A(uj) can 
alternatively be evaluated from an expansion of the en- 



ergy E at ^ = 0. 



Ox f 
Using the relation — ~ — 
ot, V ox 



it 



is straightforward to determine with the aid of (jS)) an 
expansion for the scaled energy 
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and the digamma function ^p. The zero- and first-order 
terms of the expansion l|13() are given in Using 
Eqs. (HSJ) and ^ 
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is obtained with ^lof = 



In Fig.Elthc 4th-, 5th- 



^ho.rcf 

, and 6th-ordcr expansions are compared to the results 
obtained with the non-approximated term (all for a^^ = 
—2030 ao) and with the correct atom-atom interaction 
result. Note, Eq. (|15|) can also be used for the evaluation 
of gdw, ttsc), if in the denominator S^^ei is replaced by 

p ^sc.ref 

?rof = '■ • 

^ho,ref 
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FIG. 17; Investigation of the convergence of the series expan- 
sion for the enhancement factor /^'"'"'^"(a;) in the constant 
regime. Shown are the results of the full (infinite order, dotted) 
and the 4th (dashes), 5th (chain), and 6th (dot-dash-dash) order 
expansion. For comparison, also the enhancement factor obtained 
writh the molecular potential is shown (solid). 



V. DISCUSSION AND OUTLOOK 

In this work the influence of a tight isotropic harmonic 
trap on the photoassociation process has been investi- 
gated for alkali atoms. It is found that for most of the 
states (the ones in the constant regime) there is an iden- 
tical enhancement as the trap frequency increases. This 
enhancement can reach 3 orders of magnitude for trap 
frequencies of about 100 kHz as they are reported in liter- 
ature. While the enhancement itself agrees at least qual- 
itatively with the concept of confinement of the initial- 
state wave function, also trap-induced suppressed pho- 
toassociation is possible. In fact, simple sum rule 
confirms, any enhancement must be accompanied by sup- 
pression. The physical origin of this suppression is the 
trap-induced confinement of the initial-state wave func- 
tion of relative motion within a radius that is smaller 
than the mean internuclear separation of the least bound 
vibrational states in the electronic target state. Since in 
the present calculation both initial and final state are ex- 
posed to the same harmonic trap, this result may appear 
surprising. While the explanation is based on the differ- 
ent long-range behaviors of the two involved electronic 
states, the effect itself may be very interesting in terms 
of, e.g., quantum information. 

Consider for example an optical lattice as trapping po- 
tential. The initial (unbound) atom pair is (for sufficient 
trap depths) located within a single lattice site (Mott in- 
sulator state). In the photoassociated state it could, how- 
ever, reach into and thus communicate with the neighbor 
site, if the lattice parameters are appropriately chosen. 
Such a scenario could be used for a controlled logical op- 
eration (two-qubit gate) like the CNOT. Since the latter 
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forms together with single-qubit gates a universal gate, 
this could provide a starting point for a quantum com- 
puter. Alternatively, it may be interesting to use the 
fact that if a single spot with the dimension of the trap 
length Oho or a specific site in an optical lattice can be 
addressed, then the atoms would only respond, if they 
are in their (unbound) initial state. If they are in the 
photoassociatcd excited state, they would on the other 
hand be located outside the trap and thus would not 
respond. For this it is already sufficient, if they are (pre- 
dominantly) located in the classically forbidden regime. 
Also, modifying the trap frequency it is possible to block 
the photoassociation process on demand. The trap fre- 
quency is then varied in such a fashion that a specific 
final state resonantly addressed with a laser with suffi- 
ciently small bandwidth belongs either to the constant 
or to the cut-off regime. 

A further important finding of this work is that the in- 
fluence of a tight trap on the photoassociation spectra (as 
a function of the final vibrational state) for different alkali 
atoms is structurally very similar, independent whether 
photoassociation starts from the singlet or triplet ground 
state. Also the type of interaction (strong or weak as well 
as repulsive or attractive) does not lead to a substantial 
modification of the trap influence. The only exception is 
a strong repulsive interaction that leads to a pronounced 
window in the photoassociation spectrum. The reason 
is the position of the last node in the initial-state wave 
function that in this case is located at a relatively large 
value of R and leads to a cancellation effect in the overlap 
with the final state. The nodal position depends strongly 
on Osc, but only for very tight traps also on uj. As has 
been discussed previously |43. li^ . the position of the 
window may be used for a scattering-length determina- 
tion. This will also approximately work for not too tight 
traps, but the trap influence has to be considered for 
very tight ones. Alternatively, the window provides a 
control facility, since the transition to a single state can 
be selectively suppressed. In very tight traps this effect 
is not only more pronounced, but in addition the transi- 
tions to the neighbor states are further enhanced. This 
could open up a new road to control in the context of the 
presently on-going discussion of using femtosecond lasers 
for creating non-stationary wave packets in the electronic 
excited state m One of the problems encoun- 

tered in this approach is the difficulty to shape the wave 
packet, since the high-lying vibrational states that have 
a reasonable transition rate are energetically very closely 
spaced and thus the shape of the wavepacket is deter- 
mined by the Franck-Condon factors that cannot easily 
be manipulated but strongly increase as a hmction of v. 

In view of the question how to enhance photoas- 
sociation or related association schemes (like Raman- 
based ones) the investigation of the enhancement fac- 
tors g^{i-Lj, Osc), especially its value in the constant regime 
(f/c(w,asc)) is most important. It shows that not only in- 
creasing the tightness of the trap (enlarging w) leads to 
an enhancement of the photoassociation rate, but a sim- 



ilar effect can be achieved by increasing the interaction 
strength |asc|. Most interestingly, these two enhancement 
factors work practically independent of each other, i. e. 
it is possible to use both effects in a constructive fashion 
and to obtain a multiplicative overall enhancement fac- 
tor. For a 100 kHz trap and a scattering length |asc| of 
the order of 2000 an enhancement factor (uniform for all 
states in the constant regime) of 5 to 6 orders of mag- 
nitude is found compared to the case of a shallow 1 kHz 
trap and joscl = 0. 

A comparison of the results obtained for the correct 
atom-atom interaction potential with the ones obtained 
using the approximate pseudopotential approximation or 
ignoring the interaction at all shows that these approxi- 
mations yield only for the transitions to very high lying 
vibrational states a good estimate of the photoassocia- 
tion rate. Nevertheless, despite the complete failure of 
predicting the rates to low lying states, these models al- 
low to determine the enhancement factor in the constant 
regime. For weakly interacting atoms (small |asc|) al- 
ready the pure harmonic-oscillator model (ignoring the 
atomic interaction) leads to a reasonable prediction of 
the trap-induced enhancement factor fd'^)- 

It is important to stress that the results in this work 
were obtained for isotropic harmonic traps with the same 
trapping potential seen by both atoms. In this case 
ccntcr-of-mass and relative motion can be separated and 
in both coordinates an isotropic harmonic trap potential 
(with different trap lengths due to the different total and 
reduced masses) is encountered. As is discussed, e.g., 
in l47j | where a numerical and an analytical solution 
are respectively derived for the case of atoms interacting 
via a pseudopotential, a similar separation of center-of- 
mass and relative motion is possible for axially symmetric 
(cigar or pancake shaped) harmonic traps. 

In reality, the traps for alkali atoms are of course not 
strictly harmonic. Since the present work focuses, how- 
ever, for the initial atom pair on the lowest trap induced 
state the harmonic approximation should in most cases 
be well justified. Independently on the exact way the trap 
is formed (for example by a far off-resonant focused Gaus- 
sian laser beam or by an optical lattice) , the lowest trap- 
induced state agrees usually well with the one obtained 
in the harmonic approximation, if the zero-point energy 
is sufficiently small. This requirement sets of course an 
upper scale to the applicability of the harmonic approx- 
imation with respect to the trap frequency. If lo is too 
large, the atom pair sees the anharmonic part of the trap. 
(Clearly, the trap potential must also be sufficiently deep 
to support trap-induced bound states, i. e. to allow for 
Mott insulator states in the case of an optical lattice). 

An additional problem arises from the anharmonicity 
of a real trap: the anharmonic terms lead to a non- 
separability of the relative and the center-of-mass mo- 
tion. In fact, a recent work discusses the possibility of 
using this coupling of the two motions for the creation of 
molecules 0|. Again, a tighter trap is expected to lead 
to a stronger coupling and thus finally to a breakdown of 
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the applicability of the harmonic model. 

For the final state of the considered photoassociation 
process there exists on the first glance an even more se- 
vere complication. Usually, the two atoms will not feel 
the same trapping potential, since they populate different 
electronic states. In the case of traps whose action is re- 
lated to the induced dipole moment (which is the case 
for optical potentials generated with the aid of lasers 
that are detuned from an atomic transition), the two 
atoms (in the case of Li the ones in the 2^S and the 
2^P state) will in fact see potentials with opposite sign. 
If the laser traps the ground-state atoms, it repels the 
excited ones. In the alternative case of an extremely 
far-off resonant trap the trapping potential is propor- 
tional to the dynamic polarizability of the atoms. In 
the long- wavelength limit (as is realized, c. g., in focused 
CO2 lasers |53) the dynamic polarizability approaches 
the static one, limA^oo Q;(A) = ast- The static polariz- 
abilities do not necessarily have opposite signs for the 
ground and the excited electronic state of an alkali atom, 
but in many cases different values. Then the trapping 
potentials for the initial and final states of the photoas- 
sociation process are different. The Li system appears to 
be a counter example, since for ^Li2 the average polariz- 
ability of the a'^E+ (2s-|-2s) state is predicted to be equal 
to S = azz = Oixx = 2ao(2s) = 2 x 165 = 330 ao For the 
I'^Eg {2s+2pz) state one has azz = ao(2s) -t- azz{2pz) = 
285 ao and a^x = ao(2s) -I- axx{2pz) = 292 ao yielding an 
average polarizability a « 290 ao !l3 • Thus the trapping 
potentials are expected to be very similar. This is not the 
case for, e.g., *^Rb2 where the average polarization for 
the a^E+ state is 670 ao and for l^Ej it is 1698 ao [H. 

It was checked that the use of very different values of uj 
for determining the initial and final state wave functions 
does not influence the basic findings of the present work. 
The reason is simple. Besides the very least bound states 
(and of course the trap-induced ones) the final states are 
effectively protected by the long-range interatomic po- 
tential from seeing the trap. However, if the two atoms 
are exposed to different trap potentials, a separation of 
center-of-mass and relative motion is again not possible, 
even in the fully harmonic case (a fact that was, e.g., 
overlooked in [lal)- One would again expect that this 
coupling increases with the difference in the trap poten- 
tials of the two involved states. A detailed study of the 
consequences of the coupling of center-of-mass and rela- 
tive motion due to various reasons like an anharmonicity 
of the trap or different trapping potentials seen by the 
involved atoms is presently underway. This involves also 
the case of the formation of heteronuclear alkali dimers 
where besides the occurrence of this coupling for the ini- 
tial state also a different long-range behavior of the in- 
teratomic potential has to be considered. 

Different interaction strengths occur naturally for dif- 



ferent alkali atoms as is well known and also evident from 
the explicit examples of ^Li, ''Li, and '^^K that were dis- 
cussed in this work. According to the findings of this 
work the choice of a proper atom pair (with large |asc|) 
enhances the achievable photoassociation yield quite dra- 
matically. Clearly, for practical reasons it is usually 
not easy to change in an existing experiment the atomic 
species, since the trap and cooling lasers are adapted to 
a specific one. In addition, the naturally existing alkali 
species provide only a fixed and limited number of inter- 
action strengths. 

The tunability of the interaction strength on the basis 
of Feshbach resonances, especially magnetic ones, marked 
a very important corner-stone in the research area of ul- 
tracold atomic gases. The findings of the present work 
strongly suggest that this tunability could be used to 
improve the efficiency of photoassociation (and related) 
schemes. However, it has to be emphasized that it 
is not at all self-evident that the independence of the 
scattering-length variation and the one of the trap fre- 
quency as it occurs in the model used in this work is 
applicable to (magnetic) Feshbach resonances. Further- 
more, the present work considered only the single-channel 
case while the proper description of a magnetic Feshbach 
resonance requires a multi-channel treatment. Notewor- 
thy, a strong enhancement of the photoassociation rate 
by at least 2 orders of magnitude while scanning over a 
magnetic Feshbach resonance was predicted on the basis 
of a multichannel calculation for a specific ®^Rb reso- 
nance already in |55l | . An experimental confirmation fol- 
lowed very shortly thereafter (5^. The explanation for 
the enhancement given in |55j is, however, based on an 
increased admixture of bound-state contribution to the 
initial continuum state in the vicinity of the resonance. 
This is evidently different from the reason for the en- 
hancement due to large values of \asc\ discussed in the 
present work. 

An extension within the multichannel formalism that 
allows for a full treatment of magnetic Feshbach reso- 
nances is presently underway. 
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